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A detailed study of the spin dynamics of the S-state ions Gd 3+ and Eu 2+ in the filled skutterudites 
Lai_ I i? I Pt4Gei2 (R = Gd, Eu) is reported. The spin dynamics is investigated directly by means of 
Gd 3+ and Eu 2+ electron spin resonance (ESR), performed at X-band (~ 9.4 GHz) and Q-band (« 34 
GHz) frequencies in the temperature intervals 8 < T < 300 K and 1.5 < T < 300 K respectively. 
The ESR parameters provide direct evidence for the vibrational behavior of the Gd 3+ ions but not 
for the Eu 2+ ions. These results are interpreted in the light of recent discussions about the spin 
relaxation in cage systems. In particular, the Gd 3+ spin relaxation in Lao.gGdo.iPt4Gei2 provide 
evidence for the existence of an extra phonon mode with an Einstein temperature 8e ~ 24 K in this 
system. The work suggests that the so-called "rattling" modes have a general and important role 
that should be taken into account in the study of spin dynamics in cage systems. 



I. INTRODUCTION 

Several families of intermetallic compounds are today 
known to crystallize in cage-like structures. A topic of 
growing interest in this broad area of research is the role 
of the lattice dynamics in the physical properties found 
in these materials. Among these cage systems, the fam- 
ily of the filled skutterudites has been attracting special 
attention triggered by the wide range of physical phe- 
nomena found in these compounds. These include exotic 
strongly correlated ground states [T] as well as a relatively 
large Seebeck coefficient, which opens the perspective of 
using skutterudites in the construction of thermoelectric 
devices [2J. 

Filled skutterudites have the general formula R a T 4 Xi2 
and crystallize in the LaFe 4 Pi2 structure with space 
group im3 and local point symmetry T/j for the R ions 
[3.. The R element is usually referred as guest, or filler, 
and resides in the large void in the framework of the 
[T4X12] host structure. Systematic analysis of the chem- 
ical bonding in skutterudite compounds [4] [5] revealed a 
fundamental role of the R ions in the stabilization of the 
crystal structure. 

It is generally accepted that in many of these com- 
pounds, a substantial change in the collective vibrational 
dynamics of the system is implied by the dynamical be- 
havior of the R ions (6j [7] . In the simplest approxima- 
tion, the R-ion dynamics is described in terms of a local- 
ized and isolated phonon mode, usually called "rattling" 
mode, and it can be fully described by a single parameter 
6e, the Einstein temperature [8l410|. However, the nature 
and role of the total vibrational dynamics has been a mat- 
ter of intense debate in the field, and there are a grow- 
ing number of experimental [7J [TTHH] and theoretical 
[T51 - H7] investigations demonstrating the inadequacy, and 
exploring possible extensions, of this scenario. Of special 
interest for the present work are current studies of the 
spin dynamics in skutterudites, and also in other cage 
systems, by means of electron spin resonance (ESR) and 
also nuclear magnetic resonance (NMR) [TBH25] , 

A recent theoretical work by Dahm et al. [26] investi- 
gated the NMR relaxation from rattling modes, includ- 



ing anharmonic ones. It was found that, even in metallic 
skutterudites, the rattling behavior may give rise to an 
important contribution to the spin relaxation by means 
of a Raman process (two-phonon coupling) which would 
coexist with the usual Korringa process (conduction- 
electron scattering). Recent NMR experiments in the 
skutterudites LaOs 4 Sbi2 [H| and LaPt 4 Gei2 [23] have 
successfully identified an additional contribution to the 
Korringa relaxation, that was analyzed in the framework 
of the above cited theory [2B] . The relaxation at the La 
site was shown to be field independent, which was inter- 
preted as a piece of evidence for a non-magnetic origin of 
this additional relaxation. 

In a recent work, the spin dynamics of Eu 2+ in the 
EuT 4 Sbi2 (T = Fe, Ru, Os) skutterudites was investi- 
gated by ESR [22] and clear signatures of the Eu 2+ dy- 
namical behavior were observed. It was found that the 
ESR linewidth AH, as a function of temperature, peaks 
at about 6 E and increases linearly at higher tempera- 
tures. Nevertheless, in contrast with the NMR experi- 
ments, the ESR results were shown to be field dependent. 
Thus, the results were ascribed to inhomogeneities of the 
spin-spin interaction, implied by the rattling modes. It 
was suggested that these inhomogeneities are quenched 
at high fields, causing the above features in AH to dis- 
appear. 

A fundamental issue in the application of ESR to probe 
the interplay of localized spins and the guest-ion dynam- 
ics in metallic skutterudites is the bottleneck effect |27| . 
The bottleneck effect in metals occurs when the con- 
duction electron (ce)-lattice (L) relaxation time (1/T C6 l) 
is slow in comparison with the conduction electron- 
localized spin (S) relaxation time (1/T ce s). Hence, after 
receiving energy from the localized spins, the conduction 
electrons instead of dissipating energy to the lattice, give 
this energy back to the localized spins. As a consequence, 
the relaxation of the localized spins is now modulated by 
the slow 1/T co l. In turn, one should always be careful in 
drawing any conclusions about a relation between effects 
in AH and effects in the exchange interaction between 
the local and itinerant spins. 

Aiming to discuss the applicability of the ESR tech- 
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nique in the study of the spin dynamics in cage systems, 
here, we present at two microwave frequencies (X- and 
Q-band) the ESR spectra of Gd 3+ and Eu 2+ (configura- 
tion J = S = 7/2, L = 0) in La,x- x R x PUGe u (R = Gd, 
Eu) for various values of x. Choosing S-state ions, we 
rule out non-trivial crystal-field (CF) effects implied by 
the local symmetry of the "filler" site [28 . In addi- 
tion, we address the difference between the substitution 
of La 3+ by small amounts of Gd 3+ (isoelectronic substi- 
tution) and Eu 2+ (hole doping). 

The filled skutterudite LaPt4Gei2 is a metallic com- 
pound which undergoes a superconducting transition at 
T sc = 8.3 K [23]. Electronic structure calculations show 
that the Fermi surface of this compound is mainly com- 
posed by Ap orbitals from the Ge atoms, with rather small 
contributions derived from the Pt 5d orbitals. Magneti- 
zation measurements show no signs for magnetic fluctu- 
ations or itinerant moments. The superconducting state, 
however, is not yet fully understood |30[ 13 1 J and we ad- 
dress the partial suppression of the superconducting state 
due to the inclusion of Gd 3+ moments. 

The paper is organized as follows: first, we present 
and discuss our experiments with the Gd doped samples. 
Then, we contrast these results with the experiments on 
the Eu-doped sample and also with EuPt4Gei2. After- 
ward, we explore our results in the context of relaxation 
from rattling modes [26 . 



II. EXPERIMENT 

Single-phase polycrystalline samples of 
La^Gd^^Ge^ (x = 0.005, 0.01, 0.05 and 0.1) 
and Lai_ a; Eu 2 ;Pt4Gei2 (x = 0.01, 0.05, 1) were prepared 
as described elsewhere |5J To verify the content 

of the S-state ions, the temperature dependence of the 
DC susceptibility of each sample was measured and 
the experimental results were fitted to a Curie- Weiss 
behavior, assuming the full effective moment /x e ff = 7.93 
Hb for both Gd 3+ and Eu 2+ ions. These measurements 
were performed in a commercial SQUID magnetometer 
(Quantum Design). The ESR measurements were 
performed in a Bruker Elexsys 500 spectrometer at both 
X-band (v = 9.4 GHz) and Q-band (v = 34.4 GHz) 
frequencies, in the temperatures intervals 8 < T < 300 
K and 1.5 < T < 300 K respectively. The Q-band 
measurements at relatively lower temperatures were 
possible because the typical applied magnetic fields close 
to the resonance field H res at Q-band (H res w 1.2 T, see 
below) used in the ESR measurements are larger than 
the upper critical field H c2 ( H c2 « 1 T at T w 7 K ) 
measured for LaPt4Ge 12 [5]. 

The samples where powdered and sieved (gritting of 
40/im) for grain size homogeneity. The size of these 
grains is still much larger than the skin depth that can 
be estimated from the resistivity measurements in the 
temperature interval of the experiment. Hence, we ob- 
served an asymmetric resonance line that was fitted as 



described in Ref. (35]. The parameters included in the 
fitting are the ESR linewidth AH, the resonance field 
H res , which defines the ESR g- values by applying the 
resonance condition hv = gfiBH res , the resonance am- 
plitude and an a = D/A parameter expressing the ratio 
between the dispersion (D) and absorption (A) of the 
microwave radiation when it probes a metallic material. 



III. RESULTS AND DISCUSSION 

A. ESR on Lai a; Gd ;c Pt4Gei2 

An ESR experiment probing a small amount of para- 
magnetic ions in a non-paramagnetic host lattice is sup- 
posed to bring important information concerning not 
only the ionic impurity state but also the host. In the 
present case, Gd is an isoelectronic substitution and is 
expected to have only a local effect on the lattice param- 
eter. 

In Fig. [T] we show some selected ESR spectra taken 
at both X-band and Q-band frequencies. In the whole 
temperature interval, we observe a single Lorentzian line- 
shape and, even in the lowest temperature of the experi- 
ment, no sign of unresolved CF effects were observed in 
the powder spectra. 

Fig. [2] shows the low-temperature A-band ESR 
linewidth (AH) for four different substitutions of Gd. 
It is clear that there are no concentration effects in the 
order of magnitude of AH . In the figure, the continuous 
line is the fitting to the usual Korringa expression a + bT. 
The residual linewidth a is larger for lower concentra- 
tions. This quantity is expected to be determined by CF 
effects, dipolar interactions and spin-spin exchange nar- 
rowing. The latter effect renders the residual linewidth 
narrower for higher concentrations and seems to be dom- 
inating over the broadening effects of the dipolar fields. 
The important result, however, is that the apparent Ko- 
rringa rate b is concentration independent. 

In the context of the bottleneck effect, we recall that 
since 1 /T cc g is directly proportional to the concentration 
of the ESR probe (Gd 3+ ions in this case), an increase 
of this concentration should slow down even more the 
apparent relaxation rate. Thus, the hallmark of a bottle- 
necked system is a concentration-dependent relaxation. 
The ESR AH should evolve as in a Korringa-like relax- 
ation AH(T) = a+@-T , where /3 is a constant, which in 
first approximation does not depend on the exchange in- 
teraction, x is the concentration of the ESR probe and a 
is the residual linewidth AH(T = 0). The physics of the 
bottleneck effect was well explored in the classical ESR 
literature 33 35 and also in the review by Barnes |27| . 

Our result indicates that the system is non- 
bottlenecked even for relatively high concentrations. 
From our previous discussion, one can realize that there 
are two possible routes for opening a bottleneck, being ei- 
ther increasing 1/T co l or lowering 1/T ce s. In LaPt 4 Ge 12 ; 
the conduction electrons at the Fermi surface originate 
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Figure 1: (Color online) Selected a) X-band and b) Q-band 
ESR spectra of the derivative of the resonance absorption 
(dP/dH) for La i _ ;c Gd l; Pt4Gei2 (x = 0.1 and x = 0.05). 
A single asymmetric Lorentzian lineshape is observed in the 
whole temperature interval covered by the experiment. 



mainly from the Ge 4p electrons, which are more strongly 
coupled to the lattice than s-band electrons would be. 
This contributes to the increase of l/T ce L- On top of this 
effect, one should also consider that the modification of 
the total vibrational dynamics of the system, as implied 
by the partial substitution of La by Gd, contributes to 
the scattering of the conduction electrons thus also in- 
creasing l/T ccL . 

In the non-bottlenecked regime the spin relaxation of a 
paramagnetic ion in a metallic host should be described 
by the Korringa process. The expected value for the Ko- 
rringa rate (for a g-independent exchange coupling and 
in the absence of multi-band effects) is given by [27_ : 
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Figure 2: (Color online) Low-temperature X-Band ESR 
linewidth (AH) for different substitutions of Gd in 
Lai_ :E Gd :E Pt4Gei2. The thick line is the linear fitting to the 
expression AH(T) — a + bT and the obtained parameters are 
shown in the figures. 



the full expression for both quantities |27j : 

Ag = ( V (E F )J(q = 0)) Av = V (E F ) (J(q = 0)) Av 



b=^(((N(E F )J(k F ,k' F )) 2 ) ) 

g[lB \ I Av 



(2) 



^ V (E F f { J{ q f) 



Av 



-Aff 2 



(1) 



where Ag is the so-called (/-shift, which occurs in met- 
als due to the exchange interaction between the localized 
spins and the conduction electrons [27] . A g-shift is cal- 
culated in reference to the g- values of experiments in in- 
sulators (Ag = <7 e xp-<7ins)- Adopting the low-T value of 
the g-shift, we have Ag = 0.009(4). The obtained value 
for the Korringa rate (b rj 2.3 Oe/K) is much lower than 
the one observed. Going through the relation between b, 
Ag and the exchange interaction in more detail, we write 



irki 
9Hb 



v(e f ) 2 JI 



(3) 



where in equations [2] and [3] the brackets denote an av- 
erage over the Fermi surface and rj(E F ) is the density 
of states for a given spin direction at the Fermi sur- 
face (states eV _1 mol _1 spin _1 ). From these, one can see 
that Ag is a homogeneous polarization of the Gd 3+ spins 
due to the exchange interaction with the conduction elec- 
trons, whereas b is related to a scattering process. Thus, 
in writing Eq. [T] we have assumed J\ — J2 which does 
not seem to be the case in our experiment. 
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As extensively discussed in the literature [2"Tl [551 155] , 
this contrast, in the absence of a bottleneck, may have 
its origin in multi-band and/or electron-enhancement ef- 
fects. However, previous works on the host compound do 
not support either the presence of different electronic con- 
tributions to the Fermi surface or significant electronic 
correlations [59]. We shall also discuss this issue later 
in this section, in connection with our Q-band measure- 
ments. 

A general analysis of the bottleneck phenomenon [27] 
shows that, if the spin-orbit scattering due to the mag- 
netic impurity (meaning the dependency of 1 /T ce L on the 
amount of Gd 3+ spins) is stronger than the effective ex- 
change scattering (that favors the Korringa relaxation), 
the system is non-bottlenecked. Pursuing this line, we 
suggest that, even if the homogeneous polarization (given 
by Eq. |2| is small, the exchange scattering process could 
be enhanced by the rattling modes. In this interpreta- 
tion, the rattling modes are involved in the exchange scat- 
tering process which still dominates the relaxation. 

The suppression of the superconducting state (AT SC = 
T sc ( x — 0) — T(x)) by a given concentration (Aa;) of 
magnetic impurities is also given by an effective exchange 
scattering usually comparable with the exchange scatter- 
ing given by the ESR [35 . The process is described by 
the Gorkov-Abrikosov expression: 



Aa; 8kb 

where J g is the effective exchange scattering, gj is 
the Lande gyromagnetic factor, and J is the impurity 
total angular momentum (here J — S, since for Gd 3+ 
L = 0). In Fig. [3] we show the rather small suppression 
of Tsc by the effect of magnetic impurities. The data was 
obtained from magnetization measurements measured in 
warming in 20 Oe field after zero field cooling. Tsc was 
defined by extrapolation of the steepest slope of xCO to 
x(T) = 0. Using the calculated density of states at the 
Fermi level T)(E F ) = 13.4 states eV _1 f.u. -1 = 15.8 
states eV _1 mol _1 spin _1 in Eqs. [3] and [4] we obtain 
the following estimate for the effective exchange scatter- 
ing: J E ff SR = 1.54 ± (0.40) meV and J c s £ = 1.24 ± (0.20) 
meV, which agrees reasonably well with Jeff R and make 
a strong case in favor of a relaxation process governed by 
the exchange interaction. It is also a piece of evidence 
that the relaxation is not in a bottleneck regime. Nev- 
ertheless, we still lack a mechanism for the enhancement 
of J2 = J c ff in comparison with J\ from Eq. [2] 

Figure |4ja) shows the A-band AH for the two high- 
est concentrations (x — 0.05 and x = 0.1) which allows 
one to follow the temperature evolution of the resonance 
spectra up to higher temperatures, still with relatively 
high resolution. It is shown that AH deviates from 
the apparent linear behavior at low T. Furthermore, 
the data suggest that the high temperature relaxation 
regime has a weak dependence on the concentration of 
the ESR probe. In previous ESR experiments in metal- 
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Figure 3: (Color online) Suppression of the superconducting 
state in Lai_ :E Gd a ;Pt4Gei2 by the inclusion of Gd impurities. 
The thick line is the linear fit to our experimental data. The 
value of T = 8.27 K is the transition temperature (Tsc) for 
the pure LaPt4Gei2 compound. 

lie cage systems, d(AH)/dT was also shown to be non- 
constant [18, 20 . In particular, it was found that the de- 
viations from the linear behavior may appear as a "kink" 
in AH(T), which defines a low-temperature and a high- 
temperature Korringa-like rate [20 . 

In the above cited works |18ll2"0"] , it was speculated that 
the origin of the temperature dependence of d(AH) jdT 
is due to some reorganization of the Fermi surface im- 
plied by structural effects and its concomitant changes 
in the band structure [181 12"0"] . In addition, it is gener- 
ally accepted that by multiple filling (filled skutterudites 
with more than one element as guest) one can induce 
a larger disorder in the lattice structure of the mate- 
rial, thus causing a broadening in their phonon spectra 
[361 [57] , In this sense, the origin of the concentration 
dependence of d(AH)/dT at high temperature can be 
explained by a combination of both effects: a change 
in the disorder of the phonon spectra, due to different 
Gd concentration, would lead to a different energy scale 
where the putative reorganization of the Fermi surface 
takes place resulting in a weak dependence of the high 
temperature relaxation on the Gd concentration. How- 
ever, we would like to point out that, as implied by Eqs. 
[2] and [3] similar effects would rise due to changes in the 
exchange scattering. 

Any effect in the Fermi surface, as well in the exchange 
scattering, would also be reflected in the ESR g-values. 
Therefore, the observed change of about 1/2 in the value 
of b should also manifest itself in the g- values. In Fig. 
[3Jc) we present these data, again for the x = 0.05 and 
x = 0.1 samples. Although the result is encouraging, one 
should take it carefully, since the estimated error bars are 
about half the total variation. However, our data clearly 
suggests a downturn in the g- values. 

The vibrational dynamics of the guest ion may give rise 
to local inhomogeneities of the crystal field which could 
in principle appear in an ESR experiment as an inhomo- 
geneous broadening of the resonance [38 . This may be 
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line is a guide to the eye (see text for a discussion). 



revealed in an experiment at higher frequencies. In Fig. 
[5|a)-(b) we show our Q-band measurements. It is note- 
worthy that no such broadening takes place, demonstrat- 
ing that the resonances are homogeneous. This means 
that no static disorder is being detected in our experi- 
ment (compare the order of magnitude of AH in these 
figures and Fig. |l](a)-(b)). 

In Fig [5^a) one can note that AH broadens at low 
temperature, probably due to spin-spin interactions. Fig- 
ure [5]j a) also shows that above the temperature at which 
the low-temperature magnetic fluctuations cease to dom- 
inate AH (T > 10 K), the evolution of AH with tem- 
perature is fairly described by a linear behavior with 
b = 12 Oe/K. To highlight this difference, we present 
in Fig. ^b) the Y-band and Q-band AH for the sample 
Lao.95Gdo.o5Pt4Gei2 . In the Q-band measurements, a 
systematic variation of the g- values was not observed and 
a concentration and temperature independent g- value of 
g = 2.011(5) was found. The expected Korringarate (Eq. 
[TJ obtained from this result is b = 9 Oe/K which in turn 
indicates, within experimental error, that in the experi- 
ment at Q-band the relaxation of the Gd 3+ ions follows 
the simple picture given by Eq. [l] Therefore, our disre- 
gard of multiband effects in the analysis of the Y-band 
data is also supported by these results, since one can- 
not expect that the typical fields involved in the Q-band 
measurements would eventually suppress any multiband 
effects. 

Put together, these results favors the idea that, as far 



Figure 5: (Color online) a) Q-band ESR linewidth (AH) for 
Lai_ :E Gd :E Pt4Gei2 (x = 0.05 and x = 0.1). The thick line 
is the linear fit to the expression AH(T) — a + bT and the 
obtained parameters are a — 195(1) Oe and b — 12(1) Oe/K. 
Again, no concentration dependency of AH (apart from the 
low-temperature broadening) is observed, (b) Q-band and 
X-band data for Lao.95Gdo.osPt4Gei2. 



as ESR is concerned, the rattling behavior manifests it- 
self through a coupling with the exchange interaction. A 
more comprehensive discussion on this topic will be given 
by the end of section C. Now, we point out that it has 
been recognized |39| that the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction may give rise to inhomoge- 
neous broadening and shift of the ESR spectra. Hence, it 
is plausible that an intricate relation between the Gd 3+ 
rattling and the RKKY interaction between the Gd 3+ 
spins, explains our results. 

We mention again that such field dependence of the 
relaxation was not reported in the NMR experiments 
[191 |2"5] , We also stress that the absence of static in- 
homogeneities is an important sign of the sample qual- 
ity, in particular, it demonstrates that the Gd ions are 
homogeneously distributed in the sample. 



G 




10 20 30 40 50 



Figure 6: (Color online) Jf-band ESR linewidth (AH) for 
Lai_^Eu :E Pt4Gei2 (x = 0.01 and x = 0.05). The thick line is 
the linear fit to the expression AH(T) = a + bT. For x = 0.01 
the obtained coefficients are a — 977 Oe and b = 10.5(5) Oe/K 
whereas in the case x = 0.05 the coefficients are a — 834(5) 
Oe and b — 22.5(3) Oe/K. The inset shows a temperature 
independent gr-value of about g = 1.975(5), for x — 0.01, and 
g = 1.963(5) for x — 0.05, implying in a negative g-shift (see 
text for discussion). 



B. ESR on Lai ^Eu :E Pt4Gei2 and EuPt4Gei2 

The partial substitution of La for Eu adds one hole 
to the Fermi level per f.u.. Thus, it is expected that 
this local electronic inhomogeneity gives rise to interest- 
ing effects in the interaction of the Eu 2+ spin with the 
conduction electrons. 

We present in Fig. [6] the X-band ESR AH for 
Lai_ a; Eu a ;Pt4Gei2 (x = 0.01 and x = 0.05). Again, no 
signs of a bottleneck effect were found. In comparison 
with our previous discussion, the huge residual linewidth 
(a s» 1000 Oe) is most likely due to relatively larger CF 
effects [27 . In the inset of Fig. [6] the thick solid lines 
represent the average of the measured g-values for each 
sample. From these values, one obtains relatively large 
Ag = -0.018(5) and Ag = -0.030(5) for the x = 0.01 
and x = 0.05 samples, respectively. These ^-shifts (c.f 
Eq. [TJ give a Korringa rate of b = 9(3) Oe/K and 
b = 28(8) Oe/K, respectively. In the former case, the 
measured Korringa rate is b = 10.4(5) Oe/K and in the 
latter case b — 22.4(5) Oe/K. Therefore, a comparison of 
the results suggests that for both x = 0.01 and x = 0.05, 
the Eu 2+ relaxation fits in the simple picture of Eq. [l] 

A negative Ag originates from a covalent mixing be- 
tween localized and itinerant states |271 140] . In the case 
of 4/ impurities, a negative Ag results from the mixing of 
4/ states and d itinerant electrons. Thus, our results, for 
both the Korringa rates and <?-shifts, ask for a significant 
change of the local electron density at the Eu impurities, 
even for small substitution of La by Eu. 

Different from the Gd case, the ternary EuPt4Ge 12 
compound was available for studies. Figure [7] presents 



our results. Here we find Ag = -0.035(5) (for T > 20 K), 
implying, by Eq. [T] a Korringa rate b ~ 34 Oe/K. The 
measured Korringa rate for the concentrated EuPt4Gei 2 
system is b = 8.7(3) Oe/K. This enormous contrast is 
usually attributed to the q-dependency of the exchange 
coupling J c ff, and provides evidence for an indirect ex- 
change interaction between the localized spins of the 
RKKY type gT]. 

Alternatively, it mighty be that for the concentrated 
system (x — 1) the bottleneck regime is eventually 
reached, resulting in a slowing down of the apparent re- 
laxation rate. However, in favor of the presence of an 
RKKY interaction, we point out that the negative Ag 
found in this case (Ag = —0.035(5), x — 1) cannot be 
due to changes in the local electron density, but should re- 
flect an important distinction between the Fermi surfaces 
of LaPt4Gei2 and EuPt4Gei2. Hence, since the presence 
of d electrons favors the spin-orbit scattering (and thus 
1/IceL), the system should not be bottlenecked or, at 
least, not severely bottlenecked. As an important conse- 
quence, the evolution of AH is still due to the exchange 
interaction between the local and itinerant spins , which 
is now g-dependent. 

Figure [7] also reveals a low-temperature increase of 
AH. EuPt4Gei2 is known to order magnetically below 
1.78 K. Specific heat and dc-susceptibility display four 
anomalies which indicate a complex magnetic phase dia- 
gram [323 ■ Due to the sensitivity of the ESR technique, 
it is possible that the broadening of AH is related to the 
magnetic fluctuations which give rise to these different or- 
dered states. The inset shows that the g-value increases 
for T < 20 K, most likely due to the onset of the internal 
fields associated with the discussed magnetic phases at 
low temperatures. 

In the whole temperature interval investigated, no in- 
dications could be found for an Eu 2+ rattling behavior. 
The results of our Q-band investigation are very similar 
and do not add to our discussion. 



C. Possible role of a Raman relaxation 

Following the discussion of Dahm et al. [26 , it is 
tempting to look for other sources for the relaxation 
of the localized ions, specially in the case of the Gd- 
substituted samples. It is possible that the substitution 
of even a small amount of La by Gd, gives rise to a sig- 
nificant change in the guest vibrational dynamics. In 
this context it has to be mentioned that a ternary phase 
GdPt4Gei2 does not exist as a stable phase under am- 
bient conditions. Indeed, some threshold requirements 
for the filled atoms in order to stabilize a skutterudite 
have to be fulfilled. Several studies have been devoted to 
probe high-pressure synthesis routes to compounds that 
are based on cations with radii below the critical val- 
ues (see [5] and references therein). However, attempts 
to synthesize GdPt4Ge 12 by similar methods failed. In- 
stead, phases with different crystal structures form. 
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Figure 7: (Color online) X-band ESR linewidth {AH) for 
EuPt4Gei2 . The thick line is the best fit to the expression 
AH(T) = a + bT, and the obtained coefficients are a = 412(2) 
Oe and b = 8.6(1) Oe/K. The line describes the behavior 
of AH for T > 50 K. In the inset, the ESR g-values are 
presented. For T > 20 K, g = 1.958(5) is observed, implying 
in a negative g-shift. Below this temperature the g-values 
increase (see text for discussion). 

The theory of Dahm et al. was written in terms of 
NMR parameters and does not allow to make a direct 
comparison with our results. Nevertheless, given that 
the methods for the study of the relaxation processes are 
quite general, we expect that the qualitative aspects of 
the theory will remain valid. Therefore, we simply add 
a scaling parameter to match the theoretical curve with 
our results. 

This theory considers that the guest vibrational dy- 
namics may be described by a Hamiltonian for the an- 
harmonic oscillator, which is treated in a quasi-harmonic 
approximation. This treatment gives rise to a non-linear 
equation for an effective frequency which, now, is tem- 
perature dependent: 

/ Uo\ _ 1+Q ,^oo / 1 1 1 \ 

\uooJ \cxp(haj /k B T) - 1 2 2w o/ 

(5) 

where ljq is the effective phonon frequency, a > is 
a dimcnsionless parameter characterizing the amount of 
anharmonicity and woo = w at T = 0. In the sim- 
plest rattling approximation, ujq = loqq = 0e- The relax- 
ation rate \jT\T which emerges from this picture gives 
a fair description of the relaxation phenomena in some 
skutterudites |19l |2"5] . It puts forward the idea that the 
rattling modes have a fundamental role in the spin re- 
laxation of cage systems by inducing an indirect phonon 
relaxation, usually found in experiments with insulators. 

In Fig. [8] we compare our experimental results with 
the above presented theory. In the NMR experiments 
[Ti?l |2"3"] , the contribution ascribed to the rattling modes 



was isolated after a comparison between the relaxation 
rates of the nuclear spins of the elements in the cage 
structure, and that of the rattler element. However, since 
the lattice dynamics of filled skutterudites is determined 
by significant hybridization of vibrational states of the 
guest-atoms and of the host structure [6j|7], such a par- 
titioning into independent contributions is not straight- 
forward and therefore should be only a first approxima- 
tion. Here, even such an approximate "two-site" inves- 
tigation is not possible. Our choice is to subtract the 
Korringa rate obtained at Q-band from our linewidth 
data taken at Y-band. We use the results in the case 
of Lao.gGdo.iPt4Gei2, which are obtained in a broader 
temperature range. In doing that, we assume that the 
rattling effects in the relaxation are quenched at Q-band 
and that at Y-band one may find both contributions. 
The value of woo = 24 K found in this way is consid- 
erably lower than 9e = 96 K obtained from the specific 
heat of LaPt4Gei2 [55] • This result should be interpreted 
as an extra Einstein mode associated with the Gd 3+ site, 
since it is exclusively the Gd 3+ relaxation which is being 
probed. It is also expected that, since Gd 3+ is heav- 
ier and significantly smaller than La, it would vibrate 
with a lower frequency. We also used a = 5.2, which is 
larger than previously reported in the NMR experiments 
[191 12"5] . reflecting a stronger anharmonicity due to the 
inclusion of the Gd 3+ ions in the oversized cage. We 
can speculate that the effects of a fully harmonic vibra- 
tional dynamics are averaged out and cannot be observed 
in an ESR experiment, as in the previous section where 
we addressed the results for the Eu 2+ ESR. Recent heat 
capacity measurements, taken at high fields to suppress 
the superconducting phase, reveal a phonon contribution 
compatible with an Einstein phonon with 0e ~ 20 K, in 
good agreement with our results from ESR [43] . 

It is important to make an attempt to reconcile the 
analysis of sections A and C. Dahm et al [26J considered 
that fluctuations in the electron field gradient implied by 
the rattling behavior provide the mechanism coupling the 
rattling modes to the relaxation process. In the present 
case, it is known that the ESR relaxation of a system 
of spins coupled by an interaction of the RKKY type re- 
quires the evaluation of the dynamical part of the RKKY 
interaction [39]. In this molecular field approxi- 

mation for the exchange coupling among the spin of the 
local system should be applied. This molecular field takes 
into account the finite lattice spacing between the Gd 
atoms and is determined by the lattice average of the 
RKKY interaction. Hence, the relaxation of the spin 
system will be quite sensitive to the lattice parameters. 
In this direction, we can formulate a scenario where the 
change of the total vibrational dynamics of the system, 
as implied by the dynamical behavior (rattling) of the 
Gd 3+ ions, generate local fluctuations of the RKKY in- 
teraction that are quenched at higher fields. It explains 
why the evolution of AH can be explored with a theory 
based on a Raman process and is, nevertheless, field de- 
pendent. A similar explanation should hold for the ESR 
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Figure 8: Y-band (AH(T)-a-bT) /T for La . 9 Gd . i Pt 4 Ge l2 , 
where a and b are the parameters determined in the Q-band 
measurements. The resulting data are compared with the 
solid line which originates from a theory for spin relaxation 
due to anharmonic phonons |26| . Parameters that we used 
(see text) in the calculation are shown in the figure. 



on EuM 4 Sbi2 (M=Fe, Ru, Os), where the evolution of 
AH was also shown to be field dependent in a similar 
way. 



IV. CONCLUSIONS 

In Lai_ x Gd :E Pt4Gei2, the relaxation process is non- 
bottlenecked even for relatively high Gd concentrations. 
As a consequence, the relaxation process should be de- 
termined by the exchange coupling between the local and 
itinerant spin systems. Besides this result, the effective 
exchange scattering estimated from the ESR measure- 
ments (a microscopic probe), agrees well with the one 
estimated from the Abrikosov-Gorkov expression for the 
suppression of T sc (as determined from the total mag- 



netization, a macroscopic probe), due to the addition of 
magnetic ions. 

In addition, in A-band it was found that the order 
of magnitude of the measured g-shifts implies a much 
smaller relaxation rate than the one observed and it was 
proposed that the guest vibrational dynamics could give 
rise to either an enhanced exchange scattering rate or to 
an additional relaxation process. 

The field dependence of the Gd 3+ ESR relaxation sug- 
gests that fluctuations of the RKKY local field, implied 
by the guest vibrational dynamics, should play a funda- 
mental role in the ESR response of these systems. These 
were shown to be small fluctuations that are quenched in 
going from A-band to Q-band measurements, at which 
one observes a relaxation process as expected in the sim- 
plest Korringa picture. An extra phonon mode, to be as- 
sociated with the vibrational dynamics of the Gd 3+ ions, 
was isolated in the Lao.gGdn.iPt4Gei2 with 9b ~ 24 K. 

For the Eu-substituted samples, no signatures for rat- 
tling behavior were found. The local lattice structure of 
the Gd 3+ ions in LaPt4Gei2 should be more disordered, 
due to significant difference in the atoms radii, than that 
of the Eu 2+ ions. In this sense, this work provides a 
microscopic insight on how the local vibrational spectra 
may change upon small replacements in the guest site. 

The evolution of AH with T in Lai_ ;c Eu x Pt4Gei2 is 
well explained within the simplest approximation of a 
single-band metal and g-independent J e g for small con- 
tents of Eu. The presence of the d electrons at the Fermi 
surface is suggested by the measured negative g-shift. 
With increasing content of Eu, larger shifts of the res- 
onance were observed, pointing to a continuous change 
of the Fermi surface upon Eu substitution. In the case 
of EuPt4Gei2, evidence supporting the Eu 2+ ions being 
coupled by the RKKY interaction was found. 

In particular, our work adds another piece of evidence 
for the role of the guest-ion vibrational dynamics in the 
spin relaxation in metallic compounds, where the Kor- 
ringa process is expected to dominate. 
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